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A.  STATEMENT  OF  PROBLEM  STUDIED 

■  The  objective  of  this  project  is  to  solve  Partial  Differential  Equations  associated  with 
descriptions  of  three-dimensional  objects  and  interactive  visual  design  for  manufacturing, 
design  and  simulation,  and  to  investigate  corresponding  computer  vision  and  image 
processing  problems. 


B.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

The  major  achievements  of  this  project  have  been  summarized  in  the  Interim  Progress  Reports.  The 
results  of  most  significance  fall  in  the  following  major  categories. 

■  Quasilinear  PDEs ,  the  Eikonal  Equation,  and  the  Image  Irradiance  Equation. 

■  Photometric  Stereo  and  Interactive  Visual  Design 

■  Discrete  Curvature 

■  Lattice  Boltzmann  Simulation  of  Flow  by  Mean  Curvature  for  computer  vision  and  image 
processing 

■  Electromagnetic  Scattering  from  Non-Smooth  Domains  and  new  Piecewise  Polynomial 
Wavelet  Bases 

For  the  Eikonal  Equation  and  related  first-order  nonlinear  equations  we  have  obtained  a  detailed 
understanding  of  solutions  near  singular  points.  The  main  results  are  that  there  is  a  strong  lack  of 
uniqueness  of  solutions  near  such  points  and  that  solutions  can  be  less  regular  than  both  the  function 
H  and  the  initial  data  of  the  problem,  but  that  this  loss  of  regularity  only  occurs  along  a  pair  of  curves 
through  the  singular  point. 

For  photometric  stereo  we  have  focused  on  reducing  the  number  of  photometric  images  required  and 
on  improving  the  surface  integration  from  the  normals.  A  novel  method  for  surface  integration  has 
been  developed.  It  is  different  then  traditional  methods,  in  that  it  lets  the  geometry  drive  the 
integration  path.  This  method  divides  the  surface  into  small  regions.  The  surface  is  integrated  over 
each  small  path. 

The  Four  Vertex  Theorem  is  a  classical  theorem  involving  the  extreme  points  of  Euclidean  curvature. 
This  theorem  was  extended  by  B.  Dahlberg  (see  Interim  Progress  report  Year  1)  to  cover  simple 
polygons  and  the  appropriate  notion  of  discrete  Euclidean  curvature.  We  have  generalized  this  to  the 
affine  setting.  The  correct  statement  involves  a  new  (discrete)  affine  curvature.  We  obtain  a  discrete 
six-vertex  theorem  (generalizing  the  well-known,  classical  and  nondiscrete  setting). 


We  have  studied  lattice  Boltzmann  models  for  advanced  image  processing  algorithms.  The  principal 
advantage  of  the  lattice  Boltzmann  method  over  traditional  techniques  for  obtaining  solutions  of  PDEs 
is  that  they  are  naturally  well  suited  for  fully  parallel  machines,  resulting  in  very  fast  codes.  Below  we 
describe  several  different  applications  of  the  lattice  Boltzmann  model  to  sonar  and  image  processing. 

In  recent  years  nonlinear  diffusion  has  become  a  powerful  tool  for  intra-region  smoothing  of  images. 
The  results  of  nonlinear  diffusion  can  be  used  to  obtain  an  enhanced  image  or  as  a  precursor  to 
higher-level  processing  such  as  image  segmentation,  shape  description,  and  object  detection  and 
tracking.  We  have  developed  a  lattice  Boltzmann  model  for  image  filtering  by  a  reaction-diffusion 
equation.  We  have  tested  our  methods  on  various  kinds  of  noisy  images  including  infrared  airborne 
radar  images  and  sonar  images.  The  results  are  promising. 

We  have  developed  a  lattice  Boltzmann  model  for  mean  curvature  motion.  We  have  overcome  the 
difficulty  of  finding  a  proper  collision  operator  in  the  lattice  Boltzmann  equation  by  using  the  diffusion 
generated  motion  by  mean  curvature  approach. 

Based  on  diffusion  generated  motion,  we  have  developed  a  shape  recovery  system.  Our  system  is 
capable  of  recovering  shapes  of  multiple  objects,  and  it  automatically  changes  topology  as  it 
searches  for  the  object  boundary. 

The  boundary  value  problem  for  the  time-harmonic  Maxwell  equations  for  a  domain  in  3D  or  2D,  can 
be  reduced  to  solving  a  layer  potential  integral  equation  on  the  boundary  (a  surface  in  3D  or  a  curve 
in  2D)  of  the  domain.  To  study  the  solutions  of  these  types  of  equations,  we  have  developed  a  new 
class  of  piecewise  polynomial  wavelet  bases. 
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